Relaxation Dynamics of Disordered Spin Chains: 
Localization and the Existence of a Stationary State 



Simone Ziraldo 1,2 , Alessandro Silva 3 , and Giuseppe E. Santoro 1 ' 2 ' 3 
1 SISSA, Via Bonomea 265, 1-34136 Trieste, Italy 
2 CNR-IOM Democritos National Simulation Center, Via Bonomea 265, 1-34136 Trieste, Italy 
3 International Centre for Theoretical Physics (ICTP), P.O.Box 586, 1-34014 Trieste, Italy 

(Dated: December 21, 2012) 

We study the unitary relaxation dynamics of disordered spin chains following a sudden quench 
of the Hamiltonian. We give analytical arguments, corroborated by specific numerical examples, 
to show that the existence of a stationary state depends crucially on the spectral and localization 
properties of the final Hamiltonian, and not on the initial state. We test these ideas on integrable 
one-dimensional models of the Ising or XY class, but argue more generally on their validity for more 
complex (nonintegrable) models. 
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Ergodicity is a fundamental concept of classical me- 
chanics: the properties of dynamical trajectories in phase 
space determine the long-time dynamics of a system and 
the description of the eventual stationary state in terms 
of statistical mechanics. The extension of these ideas 
to the quantum realm, pioneered in 1929 by J. von Neu- 
mann [HE], has motivated a great deal of recent research, 
mainly spurred by the experimental possibility of study- 
ing the nonequilibrium dynamics of thermally isolated 
quantum systems - most notably cold atomic gases in 
optical lattices (3J H] . A highly debated issue in the re- 
cent literature is the characterization of the long-time 
dynamics of a quantum system taken out of equilibrium 
by a sudden change of one of its parameters (a quantum 
quench). If an extensive amount of energy is suddenly 
injected in the system, will the resulting dynamics tend 
always to a well defined stationary state? And what is 
the statistical ensemble describing it? 

The stationary state existence, has been investigated 
both in generic systems [T] and in Hubbard-type mod- 
els [SH7]. A fast dynamical relaxation was recently 
observed experimentally 8 in a system of cold atoms 
and its long-time stationary state results were compati- 
ble with the generalized Gibbs ensemble (GGE) [91 ITU], 
where a set of macroscopic constants of motion are con- 
strained by the initial state. Once a stationary state is 
established, the integrability or nonintegrability of the 
dynamics appears to be the crucial ingredient: while in- 
tegrable systems in the thermodynamic limit are often 
described by a GGE [Tl"lfT3] . it is generally expected 
that the breaking of integrability will lead to thermal- 
ization [13] . 

While the relaxation dynamics of uniform systems is 
well understood, recent studies hinted towards nontriv- 
ial effects due to the breaking of translational invariance. 
In the debate two features emerged: the importance of 
distinguishing between thermodynamic limit |12[ 1 1 3] and 
finite-size effects [El US] , and the possible role played by 
localization |17j . For example, while breaking transla- 



tional invariance in the initial state could introduce cor- 
relations among different constants of motion, relevant 
for finite-size systems [El [Ej , their effect has been ar- 
gued to be negligible in predicting the stationary state 
attained by local observables [13] • While in the thermo- 
dynamic limit the breaking of translational invariance 
may not have a significant effect, localization could in 
turn play an important role, to the extent of resulting 
in the absence of thermalization even in nonintegrable 
spin chains [17] . This observation appears to be consis- 
tent with earlier numerical analysis in disordered Ising 
or XY spin chains (characterized by localization of the 
eigenstates) , where a discrepancy between the expected 
GGE and the effective stationary state was observed |16j . 

The purpose of this work is to characterize the long- 
time dynamics of disordered systems focusing on the ex- 
istence of a stationary state for local observables. We 
do so by studying the time fluctuations of local observ- 
ables in disordered spin chains which can be mapped to 
free fermionic systems. On the basis of analytical cal- 
culations, corroborated by numerics, we argue that the 
existence of a stationary state depends crucially on the 
spectral properties of the final Hamiltonian, and not on 
the initial state. While for final Hamiltonians with a con- 
tinuous spectrum associated to delocalized states a well 
defined stationary state is attained (in the thermody- 
namic limit), in systems with localized states (possessing 
a pure-point local spectrum) time-fluctuations of local 
observables persist. This will be generally associated to 
a failure of the GGE in describing long-time averages of 
many-body operators. While we have tested these ideas 
on integrable Ising or XY models, we will argue for their 
validity for more complex (nonintegrable) systems. 

Let us start by defining the problem. Consider a stan- 
dard quantum quench: an initial state |^!q), ground state 
of some Hamiltonian Hq, evolves under a different time- 
independent Hamiltonian H. Given an observable A, its 
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FIG. 1: Illustration of the two possible behaviors for the 
time-dependence of an operator A(t), showing the fluctua- 
tions around the time average, which decrease [case (a)] or 
persist at all times [case (b)]. Case (a) is a quench between 
disorder-free Ising chains (ho = 0.9, hf — 0.5), and case (b) a 
quench from a disorder-free to a disordered XX chain (with 
e = 0.3, see text for details). The inset shows an ambiguous 
case: a quench from a disordered to a clean XX chain. In all 
cases, A = c^Cj, the local fermionic density operator. 

average can be separated in two terms: 

A(t) = {^ \e iAt Ae- i6t \^ ) = A + 5A(t) , (1) 

where A = limT->.oo ^ Jq " dt A(t) is its time average and 
5A(t) is the fluctuating part. Denoting by |cv) the eigen- 
states of H with energy E a , and defining C a = (a\ipo) 
and A a i a = (a'\A\a), we get 

A(t) = Y J \ G ^\ 2 A aa + z l(E "'~ Ea)t C* a ,A a , a C a 




A) + / dne- %nt F A {n), 



where i^AJ is the so-called diagonal average [TJ [51 [TJ] 

which coincides with A (assuming no energy degenera- 
cies), and the fluctuating part 8A(t) has been recast as 
the Fourier transform of a weighted joint density of states 
Fa(P) = E a ^ a C* a ,A a , a C a 5(n-E a + E a ,). The be- 
havior of the fluctuating part SA(t) (see Fig. [I]) decay- 
ing to for large i, case (a), or remaining finite (with 
persistent oscillations), case (b), is strongly tied to the 
smoothness of Fa{Q)- If the many-body spectrum {E a }, 
in the thermodynamic limit, is a smooth continuum and 
the weights C* a ,A a i a C a make Fa(Q) still integrable, then 
SA(t) will decay to zero for large t, due to the destruc- 
tive interference induced by the strongly oscillating phase 
e -int (Ri emann _Leb eg g ue lemma). If, on the contrary, 



{E a } has an important pure-point spectrum part, i.e., 
delta functions associated to localized eigenstates, then 
one should expect persistent time fluctuations for certain 
operators. 

To exemplify the previous general arguments, we con- 
sider specific quenches for Hamiltonians of the Ising or 
XY type in one dimension: 

h = - £ (w+i + WW - E w . ( 2 ) 

3=1 3=1 

where L is the size of the chain, &j (/i = x, y, z) are 
Pauli matrices at site j with <r^ +1 = o - ^ . JJ, Jj and hj 
are nearest-neighbor spin couplings and transverse mag- 
netic fields. By a Jordan- Wigner transformation [TS] H 
can be rewritten as a quadratic form of fermionic op- 
erators Cj, and through a Bogoliubov rotation we de- 
fine a new set of fermionic operators which diagonal- 
ize H jT5J ITS] . Wick's theorem allows us to express 
A(t) in terms of sums of products of one-body Green's 
functions Gj 1 j 2 (t) = (ip(t)\<Xj x CjAip{t)) and Fj 1 j 2 (t) = 
(i){t)\c i ji c i j2 \'ip(t)}. Therefore, establishing that G jlh {t) 
and Fj 1 j 2 (t) approach a well-defined limit for large t (i.e., 
their fluctuations decay) allows us to make similar state- 
ments for a large class of operators, including spin-spin 
correlations erf a? , <7? <rf , and many others which may 
be local or nonlocal in terms of the Cj fermions. The 
case of quenches in a homogeneous Ising chain (J? = 1, 

Jj = 0), with hj — ho in Hq, and hj = hf ^ ho in H, 
is simple to illustrate: for L — > oo, we have SGj 1 j 2 (t) = 

c cos (2e k>i t), where C = 4(h - 

hf) and e^o/f = yjl + 2/i 0/ / f cos A; + h 2 ^ are the quasi- 

particle energies for Ho and H. Case (a) of Fig. [I] is an 
illustration for ji = j 2 (the local density), where fluc- 
tuations decay to as i~ 3 / 2 with oscillations [20\ This 
follows from the continuous single-particle spectrum and 
from the Riemann-Lebesgue lemma. A similar statement 
holds for ji 7^ j 2 and for Fj 1 j 2 (t). 

Now we turn to disordered Hamiltonians. To simplify 
the presentation we focus on quenches for transverse- 
field XX chains, JJ = Jj, where BCS terms are absent. 
When disorder is present, simple analytical expressions 
are lacking, while pure numerics leads to ambiguous re- 
sults. Figure [Tj for instance, shows results for SGjj(t) 
upon quenching from a clean Ho (with JJ = Jj = 1 and 

hj = 0) to a disordered H with JJ = Jj = 1 + er/j, 
hj = e£j and e = 0.3 (e sets the disorder strength and r)j, 
£j are uncorrelated uniform random numbers in [— 1, 1[). 
In this case, the fluctuations of Gjj(t) are quite clearly 
persisting at all times. The inset of Fig. [I] however, 
shows results obtained from the opposite quench (from a 
disordered H a to a clean H), and the result for SGjj(t) 
is now much more ambiguous (we prove below that fluc- 
tuations decrease to in such a case) . To make progress 
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we need a quantitative discrimination of the persistence 
of fluctuations and we introduce the mean squared fluc- 
tuations of Gj 1 j 2 (t): 

Sl h =Km^j\t\5G jlh (t)\ 2 , (3) 

which is when time fluctuations of Gj ± j 2 (t) vanish [case 
(a)], and finite if they persist [case (b)]. Physically, 

is the average fluctuation of the local density c^Cj. Let 
us denote by cj, = Ylj=i u j^j * ne fermionic operators 
diagonalizing H and with their energy. If H is trans- 
lationally invariant then /i is the momentum k and Ujk — 
e ikj l^fjj^ while if H is disordered localized eigen- 

functions. By expanding the cj's in terms of the c^'s, we 
get G jlja (t) = T,^ 2 e- i{€ ^~ e ^ )i u* jl ^u j2IX2 G^ 2 , where 
G MlM2 = (ipolc^c^lipo) . Assuming no energy degeneracy 
(i.e., e Ml = e M2 only if fi\ = ^2) 5Gj 1 j 2 (t) has the same 
expression as Gj 1 j 2 (t), except for the absence of the terms 
with Hi = /i2- The integrand in Eq. ([3]) is therefore 

\SG jlh (t)\ 2 = E E e^-^-^+^J* 

U jllJ.l U h^ U 3l^3 V 'j2^i (j ^1^2^ fj, 3 fj.4, ■ 

With the further assumption of no gap degeneracy (i.e., 
£ pi - e f>2 = £ p 3 ~ e M4 onl y if Mi = M3 and H2 = /i 4 , or 
Mi = M2 an d A*3 = /U4) |21) we arrive at the key result 

5 L- 2 = E KmJWJ 2 igwj 2 , (4) 

expressing ^j lj2 for a single realization as an 
eigenfunction-weighted sum of |G MlM2 | 2 . Disorder 
averages are performed after computing <5j lj2 , because 
we want to analyze the fluctuations of a given realization 
(averaging SGj 1 j 2 (t) would cancel such fluctuations). 
Notice that, while eigenfunction properties are buried 
in the C a and A aa i factors appearing in Fa(Q), a 
many-body theory of fluctuations starting directly from 
F A {n) is difficult [2T]. 

The nature of the eigenfunctions (localized versus ex- 
tended) plays a crucial role in Eq. Q. Regardless of dis- 
order, the |G AtlAl2 | 2 's sum to the total number of fermions 
Np in the initial state: 

El^ol4^o)| 2 =iVO. (5) 

Ml/^2 

Therefore, if the final eigenstates are extended, 
l 2 l u j2M2 1 2 ~ V-k 2 ; then S 2 iJo in scales to zero 
as Np/L 2 ~ 1/L for a system with a finite density of 
fermions. If H is clean we have to take care of the degen- 
eracy £fe = e_£ and the particle-hole symmetry, but still, 
see supplementary material (SM) , we can prove a bound 
<5| lj2 < 8Np/L 2 , indicating that fluctuations vanish for 



L —> 00. This generalizes the result of Ref. [T3] to cases 
where Hq might lead to a IG^^I 2 which has important 
nonvanishing contributions also for ki 7^ k^ (see SM). 
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FIG. 2: Average values of 5^ with j = L/2, Eqs. Q and Q , 
for quenches of XX chains when \tpo) is the ground state of 
a disordered Ho (with e = 0.3), while H is clean with hj = 
(triangles D — > C), or the opposite case (circles C — > D), 
see text for details. The averages are taken over 200 differ- 
ent disordered realizations and the error bar is the standard 
deviation of the distribution (not the error on the average) . 

Figure [2] shows the disorder average [<5|j]ai; as a func- 
tion of L m the two opposite situations described above, 
i.e., quenches from a disordered H to a clean H (D — > 
C), or vice versa (C — > D). In all cases 5j - is calculated 
from Eq. Q (with the extra terms due to degeneracies 
in the D — > C case). When is clean (D — > C data), 
[^ 2 j]aw scales to as 1/L , as expected from the bound 
discussed below Eq. [5J On the contrary, when H is dis- 
ordered (G — > D data), [<5 2 j]a« converges unambiguously 
to a nonvanishing quantity for L — > 00: time fluctuations 
survive at all times when H is disordered. This is in full 
agreement with the numerical results 22 obtained for 
the density after a quench into the localized phase of the 
Aubry- Andre model in one dimension. For smaller disor- 
der amplitude e, the situation is similar, except that the 
large-L plateau occurs for larger L, due to larger local- 
ization lengths. 

To better gauge the role of the localized eigenfunc- 
tions in making S 2 finite for L — > 00, we have ana- 
lyzed histograms of the quantities appearing in Eq. (141) . 
A histogram of |G MlP2 | 2 shows that while the average of 
|G MlM2 | 2 scales to zero as 1/L, see Eq. the distribu- 
tion of its values has large tails. To analyze these tails, 
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we work with logarithmic distributions, and define 



J 2 and Afj 1 j 2 remains 
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where A/",-, 



Sjlli^jlla | 2 | U j2M2l 



(6) 



finite when L — > oo. 

We have shown that microscopic operators have per- 
sistent fluctuations after a quench to a disordered H . 
Such fluctuations could be averaged out if one considers 
extensive operators involving sums over all sites. For in- 
stance, in a quench to a final disordered Ising chain, while 
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I u .?2mI i s a normalization constant, < 
M jlj2 < f, related to the inverse participation ratio [23] the local transverse magnetization o*(t) = 2G#(f) - f 
when j! = j 2 . Figure [3] shows P%(x) for two different has pers i s tent fluctuations, the corresponding extensive 

operator, the total transverse magnetization (per site) 
m z — L^ 1 J^j &ji has fluctuations which decrease to as 
L is increased, as we have verified. Physically, extensive 
operators effectively perform a self-averaging of the fluc- 
tuations SA(t), which then vanish in the L — > oo limit. 
Again, this agrees with the numerical results of Ref. [22) . 

Let us comment on the issue of thermalization, namely 
if an ensemble exists which is able to describe long-time 
averages of Ising or XY chains. One can prove [23] that 
for one-body operators (like Gj t j 2 (t) and Fj 1 j 2 (t)), inde- 
pendently of the size — even finite — and of the quench, 
time averages coincide with the corresponding GGE av- 
erage: this result can be traced back to the constraints 
imposed by the GGE. Thanks to the Wick's expansion, 
GGE averages reproduce the time averages for any A that 
can be expressed as a finite linear combination of powers 
of Cj and cj-, as long as the time-fluctuations of Gj 1 j 2 (t) 
and Fj 1 j 2 (t) vanish. On the contrary, when the disorder 
makes the fluctuations of Gj 1 j 2 (t) and Fj 1 j 2 (t) persistent, 
the time-averages of many-body operators might in gen- 
eral differ from their GGE value. We have verified that 
this is the case, for instance, for <7j<7j +r in quenches from 
a clean to a disordered XX chain. These results are in 
agreement with the numerical findings of Ref. [22j . 

In conclusion, we have shown that the spectral prop- 
erties and the localization or derealization of the eigen- 
states of the final Hamiltonian play an important role 
in characterizing the time fluctuations of observables. 
While the observations we made are based on the analysis 
of quenches of an integrable Hamiltonian we expect that 
the structure of Fa(FI) is what ultimately governs the 
fluctuations: integrability should not be crucial in this 
respect. The fact that, in our one-dimensional models, 
disorder localizes all eigenstates is also likely not crucial. 
For fermions hopping in a three-dimensional Anderson 
model, where there is a mobility edge separating local- 
ized and delocalized eigenstates, Eq. Q is still valid, and 
we expect that localized eigenstates contribute a finite 
quantity to £?jj 2 in the thermodynamic limit. Both these 
issues, however, clearly call for future investigations. 
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FIG. 3: Plot of P£(x), Eq. g with j = L/2 (triangles), 
and P(x), Eq. |7| (circles), for XX chain quenches from a 
clean Ho (with hj = 0) to a disordered H (with e = 0.3), for 
two values of the chain length L. [PT j (x)] a v and [P(x)] av 
(dashed and solid lines) are averages over different realizations 
of disorder, with the corresponding vertical lines indicating 
their mean values fdx x[Pj^j 2 (x)] av and J dx x[P(x)] av . 

sizes, and compares it with the unweighted distribution 
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We plot both single-instance distributions (solid points) 
as well as disorder average distributions [• • ■]„„, denoted 
by lines. The unweighted distribution P(x) is smooth, 
and self-averaging, and moves towards smaller mean val- 
ues when L increases. On the contrary, P™Ax) is more 
structured (single-instance distributions depend on the 
details of the weights |uj 1/tl \ 2 \uj 2 ^ 2 1 2 ), but its mean does 
not decrease with L, due to an eigenfunction reweight- 
ing of rare events with large values of |G> 1/12 | 2 . Physi- 
cally, this is quite transparent: similarly to what happens 
for the inverse participation ratio |23j . localized eigen- 
states are rather insensitive to the size, while extended 
states are. The fact that the mean J dx xP^ 



mains finite for L 



3102 



,x , re- 



oo for almost all realizations, is 



enough to conclude that Sj i j 2 stays finite in the dis- 
ordered H case. Indeed, Jensen's inequality implies 
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